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Abstract 

A new class of cylindrically symmetric inhomogeneous string cosmo- 
logical models is investigated. To get the deterministic solution, it has 
been assumed that the expansion (9) in the model is proportional to the 
eigen value i of the shear tensor cr' j. The physical and geometric 
aspects of the model are also discussed. 
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1 Introduction 

In recent years, there has been considerable interest in string cosmology be- 
cause cosmic strings play an important role in the study of the early universe. 
These strings arise during the phase transition after the big bang explosion 
as the temperature goes down below some critical temperature as predicted by 
grand unified theories (Zel'dovich et al., 1975; Kibble, 1976, 1980; Everett, 1981; 
Vilenkin, 1981). Moreover, the investigation of cosmic strings and their physi- 
cal processes near such strings has received wide attention because it is believed 
that cosmic strings give rise to density perturbations which lead to formation 
of galaxies (Zel'dovich, 1980; Vilenkin, 1981). These cosmic strings have stress 
energy and couple to the gravitational field. Therefore, it is interesting to study 
the gravitational effect which arises from strings by using Einstein's equations. 
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The general treatment of strings was initiated by Letelier (1979, 1983) and 
Stachel (1980). Letelier (1979) obtained the general solution of Einstein's field 
equations for a cloud of strings with spherical, plane and a particular case of 
cylindrical symmetry. Letelier (1983) also obtained massive string cosmologi- 
cal models in Bianchi type-I and Kantowski-Sachs space-times. Benerjee et al. 
(1990) have investigated an axially symmetric Bianchi type I string dust cosmo- 
logical model in presence and absence of magnetic field using a supplementary 
condition a = ajS between metric potential where a = a{t) and /? = /?(<) and 
a is constant. Exact solutions of string cosmology for Bianchi type-II, —VIq, 
-VIII and -IX space-times have been studied by Krori et al. (1990) and Wang 
(2003). Wang (2004, 2005, 2006) has investigated bulk viscous string cosmolog- 
ical models in different space-times. Bali et al. (2001, 2003, 2005, 2006, 2007) 
have obtained Bianchi typc-I, -III, -V and typc-IX string cosmological models 
in general relativity. The string cosmological models with a magnetic field are 
discussed by Chakraborty (1991), Tikekar and Patel (1992, 1994), Patel and 
Maharaj (1996). Ram and Singh (1995) obtained some new exact solution of 
string cosmology with and without a source free magnetic field for Bianchi type 
I space-time in the different basic form considered by Carminati and Mcintosh 
(1980). Singh and Singh (1999) investigated string cosmological models with 
magnetic field in the context of space-time with G3 symmetry. Singh (1995) has 
studied string cosmology with electromagnetic fields in Bianchi type-II, -VIII 
and -IX space-times. Lidscy, Wands and Copcland (2000) have reviewed aspects 
of super string cosmology with the emphasis on the cosmological implications 
of duality symmetries in the theory. Yavuz et al. (2005) have examined charged 
strange quark matter attached to the string cloud in the spherical symmet- 
ric space-time admitting one-parameter group of conformal motion. Recently 
Kaluza-Klein cosmological solutions are obtained by Yilmaz (2006) for quark 
matter attached to the string cloud in the context of general relativity. 

Cylindrically symmetric space-time play an important role in the study of 

the universe on a scale in which anisotropy and inhomogeneity are not ignored. 
Inhomogeneous cylindrically symmetric cosmological models have significant 
contribution in understanding some essential features of the universe such as 
the formation of galaxies during the early stages of their evolution. Bali and 
Tyagi (1989) and Pradhan et al. (2001, 2006) have investigated cylindrically 
symmetric inhomogeneous cosmological models in presence of electromagnetic 
field. Barrow and Kunze (1997, 1998) found a wide class of exact cylindrically 
symmetric flat and open inhomogeneous string universes. In their solutions all 
physical quantities depend on at most one space coordinate and the time. The 
case of cylindrical symmetry is natural because of the mathematical simplicity 
of the field equations whenever there exists a direction in which the pressure 
equal to energy density. 

Recently Baysal et al. (2001), Kihnc and Yavuz (1996) have investigated 
some string cosmological models in cylindrically symmetric inhomogeneous uni- 
verse. In this paper, we have revisited their solutions and obtained a new 
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class of solutions. Here, we extend our understanding of inhomogeneous string 
cosmologies by investigating the simple models of non-linear cylindrically sym- 
metric inhomogeneities outlined above. This paper is organized as follows: The 
metric and field equations are presented in Section 2. In Section 3, we deal with 
the solution of the field equations in three different cases. Finally, the results 
are discussed in Section 4. The solutions obtained in this paper are new and 
different from the other author's solutions. 



2 The Metric and Field Equations 

We consider the metric in the form 

ds^ ^ A^{dx^ ^ dt^) + B'^dy^ + C'^dz^, (1) 

where B and C are functions of x and t. The Einstein's field equations for a 
cloud of strings read as (Letelier, 1983) 

Gl^Ri~]^gl^~{pu^ ~\x,x^), (2) 

where Ui and Xi satisfy conditions 

u^Ui — —x^Xi = — 1, (3) 

and 

u'-x, = 0. (4) 

Here, p is the rest energy of the cloud of strings with massive particles attached 
to them, p = pp + X, pp being the rest energy density of particles attached to 
the strings and A the density of tension that characterizes the strings. The unit 
space-like vector x^ represents the string direction in the cloud, i.e. the direction 
of anisotropy and the unit time-like vector describes the four-velocity vector 
of the matter satisfying the following conditions 

giju'u^ = -1. (5) 

In the present scenario, the comoving coordinates are taken as 

(^0,0,0,^^ (6) 

and choose x^ parallel to x-axis so that 

^,0,0,o). (7) 

The Einstein's field equations ^ for the line-element JT]) lead to the following 
system of equations: 
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= (8) 



0^4 — 



B C ^V^ C'/^V-S Cy BC EC 

= PA\ (11) 
^1 _ -^14 Ai ( Bi Ci\ Ai ( Bi Ci\ _ 

^'^^^~'^\~B^^)'^\~B^~C)-^' ^^^^ 
where the sub indices 1 and 4 in A, B, C and elsewhere denote differentiation 
with respect to x and t, respectively. 

The velocity field is irrotational. The scalar expansion 0, shear scalar cr^, 
acceleration vector Ui and proper volume are respectively found to have the 
following expressions: 

, 1 (Ai Bi CA 



^ V ^ C , 

2 1 ji l/,2 1 ^4^4 , i?4C'4 , C4A4 

a =-a,,a =-9 - — + ^ + 



(14) 



u, = u^,ju^ = (^^,0,0, Oy (15) 

= = A^BC, (16) 

where 5 is the determinant of the metric ([T]). Using the field equations and the 
relations ([TB]) and p4|) one obtains the Raychaudhuri's equation as 

^ = u;, - - 2(7^ - ipp, (17) 

where dot denotes differentiation with respect to t and 

R,,v'i^ = ^Pp. (18) 

With the help of equations H]) - (H)), the Bianchi identity (t!^^ reduced 
two equations: 

and 

Ai fAi Bi Ci\ 

Thus due to all the three (strong, weak and dominant) energy conditions, one 
finds p > and pp > 0, together with the fact that the sign of A is unrestricted, 
it may take values positive, negative or zero as well. 



to 
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3 Solutions of the Field Equations 



As in the case of general-relativistic cosmologies, the introduction of inhomo- 
geneities into the string cosmological equations produces a considerable increase 
in mathematical difficulty: non-linear partial differential equations must now be 
solved. In practice, this means that we must proceed either by means of ap- 
proximations which render the non- linearities tractable, or we must introduce 
particular symmetries into the metric of the space-time in order to reduce the 
number of degrees of freedom which the inhomogeneities can exploit. 

Here to get a determinate solution, let us assume that expansion (9) in the 
model is proportional to the eigen value i of the shear tensor cr' j. This 
condition leads to 

A = (BCr, (21) 
where n is a constant. Equations ([9]) and lfTO|) lead to 



Bm -Bii _ C44 Cii 
~B B~ C~ 



(22) 



Using (121]) in (ini), yields 



■ ''''-2n(^ + ^)(^ + ^]=0. (23) 



B C \ B C J \ B C 

To find out deterministic solutions, we consider the following three cases: 
{i)B = f{x)g{t) and C = h{x)k{t), 

{ii)B ^ f{x)g{t) and C ^ f{x)k{t), 

{iii)B ^ f{x)g{t) and C ^ h{x)g{t). 

The two cases (i) and (ii) are discussed by Baysal et al. (2001) and the last 
case (iii) is discussed by Kilinc and Yavuz (1996). We revisit their solutions and 
obtain a new class of solutions for all these cases and discuss their consequences 
separately below in this paper. Our solutions are different from these author's 
solutions. 

3.1 Case(i): 

B = f{x)g{t) and C = h{x)k{t) 

In this case equation reduces to 

h/f (2n-l)(fc4/fc) + 2n(.94/.9) 



hi/h {2n~ l){gi/g) + 2n{ki/k) 
which leads to 

/ " h 



i\:(constant), (24) 
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and 

k^/k K - 2nK - 2n , , , , 

= 2nK + 2n-l = (^6) 

From Eqs. ([25]) and (|26|) . we obtain 

/ = ah"" (27) 

and 

fc = (28) 
where a and 5 are integrating constants. Eq. ([22]) reduces to 

544 _ ^44 ^ /ll _ ^ ^ |'2g\ 

g k f h 
where is a constant. Using the functional values of B and C in (|22|) . we obtain 

gg44 + agl = --J^9^, (30) 

which leads to 

g = f3^ cosh^{bt + to), (31) 
where p and are constants of integration and 



b= ^/a{a + l). 

Thus from Eq. ^ we get 

k = (5/3*r cosh^ {bt + to). (32) 
From Eqs. ^Ei) and ((291), we obtain 

hhn+Khl = -^h^, (33) 

iV — 1 

which leads to 

h = £T<+^ cosh"f^ (rx + xo), (34) 
where i and a;o a-re constants of integration and 



Hence from Eq. (P7)) we have 

/ = ai?KTT cosh"KTT (^rx + xq). (35) 

It is worth mentioned here that equations (|30p and (|33p are fundamental basic 
differential equations for which we have reported new solutions given by equa- 
tions dSTl) and iMl). 
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Thus, we obtain 

B = fg = Q cosh"?^ {rx + xo) cosh^ {bt + to), 
C ^ hk = RcoshT^ (rx + xq) cosh^(6t + to), 



and 
where 



A = [BCY = Af cosh"(ra; + Xq) cosh"(6i + tg), 



M = (gi?)". 

Hence the metric ([Ij takes the form 

ds^ = cosh^"(rx + a;o) cosh^"(6i + to){dx^ - dt^) + 

cosh'iTTT (rx + Xo) cosh^ (6t + to)c^2/^ + 
cosh'fTT(rx + Xq) cosh^(6t + ta)dz^. 
By using the following transformation 

rX — rx + Xq, 

Z ^ Rz 
bT = bt + to 

the metric ([5^ reduces to 

Af2cosh2"(rX) cosh2"(6T)(dX2 - dr2) + 
cosh"???! (rX) cosh^ {bT)dY^ + cosh"^ (rX) cosh^ {bT)dZ^. 



(36) 
(37) 

(38) 



(39) 



(40) 



(41) 



In this case the physical parameters, i.e. the energy density (p), the string 
tension density (A), the particle density (pp) and kinematical parameters, i.e. 
the scalar of expansion (9), shear tensor (cr), the acceleration vector (iii) and 
the proper volume {V^) for the model (f4T|) are given by 



+ FTif) tanh2"(6r) +r^(n+ jj^) ts^nh^rX) 



X = 



M2 cosh2"(rX)cosh^"(for) 
(" + (^) tanh2"(6T) - (ri + j^Tw) ta.nh^{rX) + i 



M2 cosh^"(rX)cosh^"(6T) 



(42) 



(43) 
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M2 cosh^" {rX ) cosh^" (bT) 
b{n + 1) tanh(6T) 



(44) 
(45) 



M cosh" (rX) cosh" (6T)' 

^2 ^ bHanh\bT)[ia+l)^in^-n + l)-3a] 

3(a + l)2Af2cosh2"(rX)cosh2"(6T) ' ^ ' 

ii^ = [nr tanh(rX), 0, 0, , (47) 

- ^ = cosh2"+i(rX) cosh2"+i(6r). (48) 
From equations (|45|) and (I46p . we obtain 

(a + l)2(n2 - n + 1) - 3a 



3(n+l)2(a + l)2 



(constant). (49) 



The models (|41l) represents expanding, shearing and non- rotating universe. If we 
choose the suitable values of constants K and M , we find that energy conditions 
P > 0, pp > are satisfied. Since ^ is constant throughout, hence the model does 
not approach isotropy. In this solution all physical and kinematical quantities 
depend on at most one space coordinate and the time. 

3.2 Case(ii): 

B = f{x)g{t) and C = f{x)k{t) 

In this case equation (^5]) reduces to 

(4n-l)4('^ + ^')=0. (50) 



f \9 k 
The equation ([5U|) leads to three cases: 

(a) " = ^' 
5 « 

The case (a) reduces the number of equation to four but, with five unknowns 
which requires additional assumption for a viable solution. In the case (b), the 
model turns to be a particular case to the Bianchi type-I model. Therefore we 
consider the case (c) only. 
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Using condition (c) in equation leads to 

544 ^44 fr^^ 

By using condition (c) in ([?T|) . we get 

g = e^^, k - e-^^, (52) 

where T = t + ^, to, b, and L are constants. From equations ([9]) or pO|) and 
([5^. we have 

Solving ([SB]) , we obtain 

/ = cosh2"+i(Mos + xq), (54) 

where 



V2n 

Mo - 



2n + 1 

and £ and xo are constants of integration. 

It is important to mention here that (|53p is the basic equation for which new 
solution is obtained as given by ([51)1 . 

Hence, we obtain 

B = fg = Qae^^ cosh2"+^ (Mqx + xo) (55) 

and 

C = fk = Qoe-^^ cosh^'^+\Mox + xq), (56) 
where Qo = Therefore 

A = [BCY = No cosh2"(2»+i) (Mo.T + .xq), (57) 

where A^o = Qo"- 

After suitable transformation of coordinates the metric H]) reduces to the 
form 

ds^ = iV2cosh2"'"(Moa;)(dX2-dT2) + Q2^osh'"(AfoX)(e2^^dF2 + e-2i^r^^2)^ 

(58) 

where m = 2(2ri + 1). 

In this case the physical parameters p, A, Pp and kinematical parameters 9, 
a, iii and for the model (|55|) are given by 

p = 5 \Un + l)(2n - 1) tanh^fMoX) - (2n + 3)1 , 

N^(2n + 1) cosh (MqA ) 

(59) 
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r2 

[{An + 1) tanh^(Mo^) - l] , (60) 



7V2cosh2""'(MoX) 



Pr, = \n{An + 1) tanh^(MoX) - (n + 1)1 , (61) 

7V2(2n+l)cosh2"™(MoX) L ^ ^ I o ; v , V ; 

6* = 0, (62) 

2 -^^ 

" iV2cosh2"'"(MoX)' ^^^^ 

= (nmMa tanh(Afo^), 0, 0, o) , (64) 

^3 - = QJT cosh™(2"+i) (MoX). (65) 



In this case the expansion 9, in model ([58|) . is zero. With the help of physical 
and kinematical parameters, we can determine some physical and geometric fea- 
tures of the model. All kinematical quantities are independent of T. In general, 
the model represents non-expanding, non-rotating and shearing universe. The 
acceleration vector ii is zero for n — 0, n = — i. Choosing suitable values for 
n, we find that energy conditions p > 0, pp > arc satisfied. The solutions 
identically satisfy the Bianchi identities given by p9|) and ((20|) . In this solution 
all physical and kinematical quantities depend on at most one space coordinate. 



3.3 Case(iii): 

B = f{x)g{t) and C = h{x)g{t) 

In this case equation reduces to 

The equation leads to three cases: 

(a) ^=\' 

ib) ^ = 0, 
9 

The case (a) reduces the number of equation to four but, with five unknowns 
which requires additional assumption for a viable model. In the case (b), which 
infers a constant g refers to the static solution. Therefore we consider the case 

f^.O, (6T) 
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which produces non-static and physically meaningful solution as follows. Equa- 
tion (|67p leads to 

Equation (|68p . after integrating, gives 

f = e^°^, h^e-^"^, (69) 

where X ^ x + ^ and xq, r, Lq are constants. From equations ^ or pO|) and 
we have 

which after integration gives 

g = £l''+\osh'"+\Kot + to), (71) 



where Kq = 2n^\ ^"^^ ^O' constants of integration. It is important to 

mention here that (ffP]) is the basic equation for which new solution is obtained 
as given by [7l|) . 

Thus we obtain 

fg^ De^°^ cosh2"+^ {K^t + to) (72) 

and 

C^hg = De-^"^ cosh^"^+\Kot + to), (73) 
where D = Therefore 

A = (BC)" = Pcosh2"(2"+i)(A'ot + to), (74) 

where P = D'^". 

After suitable transformation of coordinates the metric ([1]) reduces to the 
form 

ds^ = P'^ cosh'^'''' {KoT)idX^^dT^)+D^ cosh' (KoT) {e^^°^ dY^ + e'^^"^ dZ^) , 

(75) 

where s = 2{2n + 1). 

In this case the physical parameters p, X, pp and kinematical parameters 9, 
a, iii and for the model ([75]) are given by 



K^{An + l){2n + 1)^ tanh2(ii:or) - i§ 
P2 cosh2"''(ii'or) 



ii:^[(4n + 2) - (4n + l)(4n2 - 1) tanh^(i^o7')] + 

p2 cosh2"^(i^or) '■ 



(77) 
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pp T;:rzz:r2KFr;:r7^^ ' ^^^> 



K^[An{2n + l)(4r t + 1) tanh^(jroT) - (4n + 2)] - 21^ 

Kos{n+l)tSiBh{KoT) 



Pcosh"''(ifoT) 



(79) 



(80) 



Therefore 



12P2cosh2"'*(i^or) 

= (0,0,0,0) (81) 
= ^ {PD f cosh(2"+i)'^(i^oT). (82) 

0-2 4^2 - 4n + 1 , ^ 

^ = 12(n+l)2 = (8^) 

The models (|75p represents an expanding, shearing and non-rotating universe. 
We find that energy conditions p > 0, pp > are satisfied if we choose n> 
and Lq 7^ and we get physically significant string cosmology model. The 
energy density p, the string tension density A and particle density pp at all 
finite spatial location tend to constant value as T ^ 0. Since we observe that 
^ is constant throughout, hence the model does not approach isotropy. The 
solutions identically satisfy the Bianchi identities given by and (PO)) . In 
this solution all physical and kinematical quantities depend on at most one 
time coordinate. 



4 Concluding Remarks 

In the study, we have presented a new class of exact solutions of Einstein's field 
equations for inhomogeneous cylindrically symmetric space-time with string 
sources which are different from the other author's solutions. In these solu- 
tions all physical quantities depend on at most one space coordinate and the 
time. 

In case (i), the models (1411) represents an expanding, shearing and non- 
rotating universe and all physical and kinematical parameters depend on at 
most one space coordinate and the time. 

In case (ii), the model ()58p represents non-expanding, non-rotating and 
shearing universe. The solutions identically satisfy the Bianchi identities given 
by ([TO]) and ([^0]) . In this solution all physical and kinematical quantities depend 
on at most one space coordinate. 

In case (iii), the models (j75l) represents expanding, shearing and non-rotating 
universe. The solutions identically satisfy the Bianchi identities given by (I19p 
and (j20p . In this solution all physical and kinematical quantities depend on at 
most one time coordinate. 
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